Journal of Heredity

Rapid evolution allows populations to persist in environments where they would otherwise go extinct. This phenomenon, known as evolutionary rescue, is typically studied in the framework of biological evolution, yet adaptive traits can also arise and spread through cultural evolution. The present study developed a stochastic eco-evolutionary model to compare rescue probabilities through biological and cultural evolution. Transmission bias governed the rescue probability under cultural evolution by setting how readily a rare adaptive trait was copied. Conformity bias suppressed population persistence because a rare trait was the least likely to be copied. Content bias toward the adaptive trait enabled evolutionary rescue when social learning was rapid, but it typically yielded a lower rescue probability than biological evolution. Only anticonformity bias, together with a high social learning rate, exceeded the rescue probability of biological evolution by enabling the adaptive trait to be established more rapidly. These results demonstrate that transmission bias alters the demographic consequences of cultural evolution and highlight the importance of transmission processes in evolutionary rescue theory. Understanding how adaptive behaviours are socially transmitted may also improve predictions of animal population persistence and inform conservation efforts in rapidly changing environments.

High-phenolic plant genera present a major technical limitation in genomic research. Standard extraction approaches that perform reliably across diverse flora often perform poorly when applied to recalcitrant taxa, producing low DNA yield and integrity incompatible with sequencing requirements. The genus Anigozanthos (Kangaroo paws) from the family Haemodoraceae exemplifies this problem. We identified key physicochemical factors governing extraction failure in this genus and resolved them through targeted modifications to lysis chemistry and contaminant management. The resulting protocol achieved a near threefold improvement in DNA purity, substantially reducing contaminant carry over and consistently yielded high-integrity, long DNA fragments (DIN > 7) across a diverse sample set spanning cultivated and wild material across four diverse genera of Haemodoraceae. We also tested a straightforward purity assessment framework that can be implemented in any standard molecular laboratory, enabling rapid pre-submission quality assessment without the need for specialised equipment. Together these advances open a practical path to genomic characterisation of Anigozanthos that establishes a transferable model for genomic research across Australia ' s chemically complex native flora.

Three new species of sun orchid (Thelymitra) endemic to Aotearoa New Zealand are here described. These are T. palustris, T. scabrifolia and T. semaphora. The morphological distinctiveness of these three species has been acknowledged for decades; however, their taxonomic status has remained unresolved. Evidence from existing karyological data, recently generated DNA sequence data (LFY and ycf1) and morphological studies from historical and fresh collections are used here to support their formal description. Both, T. palustris and T. semaphora are restricted to wet habitats north of Auckland (North Island). Thelymitra scabrifolia inhabits mostly scrub, and it has a similar northern North Island distribution, but is has been found also in Manawat[a]whi / Three Kings Islands and historically in Otago (South Island). All three species are polyploids and are of conservation concern.

Extrafloral nectaries (EFNs) are a widespread plant defense mutualism trait and are highly convergent, appearing in hundreds of plant lineages worldwide. Here we investigate a report of possible EFNs in Erythranthe angulosa, a recently-described California wildflower. We integrate field observations, insect bioassays, an induction experiment, and microscopy to test for signatures of EFN function, finding no evidence that the distinctive axillary swellings produced by E. angulosa function as EFNs. We also uncovered two distinct morphs at the type locality of E. angulosa that diverge in the number of axillary swellings produced, as well as other shoot architecture traits such as stem thickness, leaf size, and branch number. Although the axillary swellings appear to not function as EFNs, they remain a compelling morphological variant within the yellow monkeyflowers that may perform storage or another unknown function.

Intertidal mudflats are low hydrodynamic energy environments hosting microphytobenthic communities that experience strong spatiotemporal variability in light regimes, including changes in spectral quality and light intensity that can lead to cellular photooxidative stress. To cope with these fluctuations, autotrophs exhibit diverse and highly plastic adaptations that are often species-dependent and shaped by their ecological niches. This study investigates photophysiological responses and metabolic remodeling in a diatom assemblage originating from a natural winter microphytobenthic biofilm under contrasting red and blue light intensities. To this end, photosynthetic parameters were monitored alongside changes in lipophilic metabolites, including untargeted lipids and lipophilic pigments. While few metabolites showed temporal remodeling, rapid and contrasting changes were observed within 30 minutes in response to both spectral quality and light intensity. Red light treatments induced broader remodeling of lipophilic metabolites than blue light, whereas blue light appeared to have a greater impact on photosynthetic parameters. Moreover, red light induced xanthophyll-cycle responses comparable to those observed under blue light at equivalent incident intensity. We discuss these metabolic responses in relation to diatom photoadaptive strategies, placing these findings within the intertidal environmental framework. This work further underlines the importance of understanding rapid metabolic plasticity in coping with light fluctuations, providing new insights into the photoregulatory strategies of natural microphytobenthic communities.

Animal mitochondrial rRNAs are commonly viewed as structurally reduced, yet sponge mt SSU rRNAs range from compact to highly expanded structures. Using nine conserved structural anchors, we compared 216 taxonomically resolved records from four classes and 22 orders, including 16 freshwater Spongillida and 200 marine sponges. Twelve homologous hypervariable substructures were coded as structural types, and their ordered combinations as composite types. We identified 38 structural types and 62 composite types across molecules ranging from 853 to 2,019 nt. Hexactinellida and freshwater Spongillida were each uniform for a distinct composite type but differed markedly in overall structure: hexactinellid mt SSU rRNAs were compact, whereas those of Spongillida were long and contained four to five candidate insertion regions. These results show that a conserved scaffold can accommodate extensive lineage-associated structural variation and provide a practical framework for comparing highly divergent mitochondrial rRNAs.

Heritable bacterial symbionts are pervasive in terrestrial arthropods, often imposing reproductive manipulations to promote their own spread within host populations. Co-infections are common, potentially allowing symbiont co-infectors to hitchhike through a host population. However, adverse thermal conditions can disrupt these communities, particularly when co-infectors vary in their thermal sensitivity. We used a multi-generation experiment to test whether warm (29 {degrees}C) conditions disrupted spread of heritable symbionts through uninfected populations of the spider, Mermessus fradeorum. We tested two common infection combinations: a single infection with a cytoplasmic incompatibility (CI) inducing Rickettsiella or a feminizing co-infection that included a feminizing Wolbachia, the same Rickettsiella, and up to three apparent hitchhikers (two additional Wolbachia strains and Tisiphia). We initiated replicate populations with 1/3 of one infection type and 2/3 uninfected spiders, evaluating population infection rate over 5 spider generations under different temperature regimes. Under cool (21{degrees}C) conditions, Wolbachia feminization drove co-infection to 88% and Rickettsiella CI drove single infection to 83% of host populations. Vertical transmission for all symbionts was high (97-99%) and hitchhiking symbionts also spread effectively. Under warm conditions, feminization and CI efficacy were reduced, and symbionts suffered variably reduced vertical transmission. Warm conditions ultimately destroyed the co-infecting symbiont consortium and impeded symbiont spread. On its own, though, Rickettsiella was still able to increase, despite reduced strength of CI. We hypothesize that contrasting tensions between feminizing spread of the symbiont consortium versus environmentally driven loss of function and transmission may explain observed patterns of mixed infections in field populations of this spider.

Aegilops peregrina is a wild allotetraploid wheat wild relative and an important source of genetic diversity for stress tolerance and agronomic traits. Here, we report a subgenome-resolved, chromosome-scale reference genome assembly of a drought tolerant and stem rust resistant Ae. peregrina accession PI 604178 generated using PacBio HiFi and Hi-C sequencing. The 10.13 Gb assembly contains 98.81% of sequence anchored to 14 pseudomolecules representing the seven S and seven U chromosomes, with contig and scaffold N50 values of 25.84 and 746.48 Mb, respectively. The assembly achieved a consensus quality value of 74.61, 97.83% k-mers completeness, and 99.9% BUSCO completeness. LTR Assembly Index values of 20.43 and 18.79 for the S and U subgenomes, respectively, further supported high continuity across repeat-rich regions. Repetitive elements comprise 85.93% of chromosome-anchored assembly. We annotated 59,910 high-confidence protein-coding genes, with comparable gene representation across the two subgenomes. This reference genome provides a high-quality genomic framework for comparative analyses, characterization of important loci regulating agronomic and resilience related traits, and sequence-guided exploitation of Ae. peregrina allelic diversity for wheat improvement.

Blumeria graminis f. sp. avenae (Bga), the causal agent of oat powdery mildew, is one of the most host-specialized members of the B. graminis species complex. Despite its agricultural importance, the lack of a high-quality reference genome has limited studies of host specialization, virulence evolution and comparative genomics in this pathogen. Here, we generated the first chromosome-scale genome assembly of Bga using an integrative approach combining long- and short-read sequencing, Hi-C scaffolding and transcriptome data. The Bga genome exhibits hallmark features of powdery mildew fungi, including extensive repeat content and low gene density. Comparative analyses revealed that genome expansion is primarily associated with historical transposable element proliferation rather than recent transpositional activity. Genome organization is consistent with a functionally stratified "one-speed" model, in which genes associated with pathogenicity, including predicted effectors and infection-responsive genes, are preferentially located in transposable element-rich regions characterized by reduced synteny conservation and extended intergenic spaces. In contrast, conserved genes are concentrated in compact genomic regions and maintain strong syntenic conservation across cereal-infecting formae speciales. Hi-C analyses demonstrated a highly structured chromatin architecture and revealed genome organization patterns associated with infection-related gene expression. Comparative genomic analyses indicated that host specialization in Bga is driven by localized diversification of a relatively small subset of genes rather than large-scale genome restructuring. These results provide the first high-quality genomic resource for Bga and offer new insights into the evolutionary mechanisms underlying host specialization in powdery mildew fungi.

Squid abundance has been reported to increase globally between 1970 and 2010. This increase has been hypothesized to result from two primary factors: the loss of top predators due to overfishing and rising ocean temperatures. The decline in apex predators may lead to the expansion of squid populations either through reduced predation pressure or diminished competition with juvenile predators. Concurrently, increased temperatures could enhance the somatic growth rates of squid, thereby accelerating their population growth. However, empirically disentangling the impacts of predator loss and temperature on squid biomass remains challenging, especially in a food-web context. In this study, we used a size- and trait-based model of upper trophic levels that resolves the ecosystem structure -- biomass and trophic interactions of fish and squid -- for varying depth, temperature, and secondary production, to investigate two hypotheses of the historical expansion of squid, i.e., the effects of predator depletion from fishing and rising temperatures on squid biomass. Our model reveals that intensified fishing of squid predators -- specifically large demersal fish in shelf systems and large pelagic fish in open oceans -- leads to a slight increase in squid biomass. Conversely, elevated temperatures are associated with a decline in squid biomass. This temperature-driven reduction in biomass is attributed to an increased metabolism of squids beyond the available food supply. If historic overfishing on large marine predators continues to be curtailed, we expect a corresponding reduction in global squid biomass and fisheries potential, which could be further exacerbated by rising temperatures.

Most evidence of the consequences of climate change in natural systems has focussed on shifts in mean temperature (1,2), but the effects of extreme climatic events (ECEs) remain far less understood. This is particularly true for very severe ECEs that may occur only once every few decades. Understanding the consequences of these severe events for natural populations is nonetheless critical, since their frequency is predicted to rise under current climate change (3). Here we combine a unique long-term dataset spanning almost five decades of breeding (>20,000 events) and morphological data (>120,000 observations) in adult and nestling great tits (Parus major) and blue tits (Cyanistes caeruleus) with fine-scale temperature records to examine the effects of an unprecedented heatwave in May 2026 on breeding success and morphology. Average temperature during the heatwave (22-29 May 2026) was 7.85 C above the historical record, reaching +10.5 C (+4.32 SD) at its peak (25-26 May). These record-breaking temperatures significantly reduced adult breeding success and nestling bmass relative to expectation in the absence of a heat-wave. Given the heatwave was widespread (Fig. 1A), our findings from a single, exceptionally well-studied population are likely to generalise to other species exposed to the same event, providing key evidence that severe ECEs can substantially harm wild populations.

The exponential increase in the number of genotyped animals, combined with the availability of high-density SNP chips has introduced computational challenges for routine genomic evaluations, particularly during the construction of the genomic relationship matrix. Although higher-density SNP panels can facilitate the identification of causal mutations, their use substantially increases computational requirements without a proportional gain in genomic prediction performance. To optimize computational efficiency while maintaining accuracy of genomic predictions, this study compared five SNP selection strategies (i.e., random sampling, random sampling with inclusion of informative SNPs, linkage disequilibrium (LD)-based pruning, a Shannon entropy-based machine learning approach, and [[EQUATION]]-based prioritization) to develop reduced-density panels for Nellore cattle. Using high-density (HD) genotype data comprising 437,650 SNPs from 304,782 animals (after quality control) as reference, three reduced-density panels (25K, 45K, and 65K SNPs) panels were tested across five traits (i.e., Age at first calving, Stayability, Weaning weight, Yearling weight, Muscling) with diverse genetic architectures. Genomic estimated breeding values (GEBVs) derived from these reduced panels were compared to those obtained from the HD reference panel using Pearsons correlations, under both genomic best linear unbiased prediction (GBLUP) and single-step GBLUP (ssGBLUP) methods. In the GBLUP model, prediction accuracy generally improved with increased marker density. Random selection with and without the informative SNPs consistently yielded the highest accuracies, whereas the [[EQUATION]]-based approach showed the lowest agreement with the HD reference across all densities. In contrast, ssGBLUP demonstrated strong robustness to marker reduction, producing uniformly high correlations {approx}1.00) across all SNP densities and selection strategies. These findings indicate that optimized low-density SNP panels maintain prediction accuracy comparable to HD panels, offering a cost-effective tool for large-scale genomic evaluations.

The ocean tides influence key Earth system processes at a range of spatial and temporal scales. It is known that the geometry of ocean basins is the leading controller of tidal energetics, so well-constrained palaeogeographic reconstructions and tidal properties for Earths past are imperative when investigating other Earth system processes. Here, we present a novel way to constrain both deep-time tidal model results and reconstructions, by combining palaeoecology with sedimentology. We compare new palaeo-tidal model simulations for the Cambrian period, significant for the early origin and radiation of major animal fauna, to tidal proxies. One of the most abundant soft-bodied organisms preserved during this time are cnidarian medusae (jellyfish). A total of 17 cnidarian medusae localities were obtained through the literature, which had an adequate global distribution and occurred at regular intervals throughout the period of study. In some locations there were also estimates of palaeo-tidal range. Our results show a good agreement between the simulations and proxy data. In the few locations where there is disagreement, it is proposed that the palaeogeographic reconstructions are missing details, e.g., island chains, and our results allow for the palaeogeographic reconstructions to be improved. The proxy method presented is promising and can be applied to other time-periods with different marine fossils, particularly at evolutionary and extinction periods where the marginal marine environment is of importance.

Sicyoidochytrium minutum DNA virus (SmDNAV) was isolated several years ago from a protist host of family Thraustochytriaceae of the class Labyrinthulomycetes. This virus shared little similarity to other viruses in gene content and protein sequences, albeit seemingly belonging to the phylum Nucleocytoviricota. By extensive searches in genomic and metagenomic sequence databases, we identified numerous long contigs related to the SmDNAV genome and analyzed proteins shared by these putative viruses. Phylogenetic analyses place these viruses within the class Megaviricetes, outside of all established orders, and as a sister group to the clade combining families Mamonoviridae and Manesviridae. Homologs of SmDNAV proteins were found in association (either integrated or co-sequenced) with other Labyrinthulomycetes and Rhodophyta protists from diverse marine and freshwater environments. Consequently, we propose SmDNAV as the prototype member of a new order, provisionally named Ariadnavirales, within class Megaviricetes, phylum Nucleocytoviricota. Members of Ariadnavirales have lost most of the genes encoding components of the replication and transcription systems that are otherwise conserved in nucleocytoviricots, suggestive of transition to genome replication and expression dependent on the host nucleus.

Accelerating photoprotective regulation to improve carbon assimilation is a promising strategy to increase crop productivity. Although rapid non-photochemical quenching (NPQ) relaxation has been validated as a target through metabolic engineering, it remains unclear whether conventional breeding has improved this trait. Here, we investigated whether more than a century of soybean breeding enhanced NPQ relaxation alongside light-saturated carbon assimilation and seed traits. We evaluated a historical panel of 24 soybean genotypes across vegetative and reproductive developmental stages by integrating NPQ relaxation, gas exchange parameters, xanthophyll-cycle pigment profiles, expression of key photoprotective genes (VDE, PsbS, and ZEP), seed number and seed weight. NPQ relaxation parameters were not consistently associated with genotype release year, seed number, or seed weight at either developmental stage. The only exception was the amplitude of the rapidly relaxing NPQ component (AqE), which was negatively correlated with all three variables during the reproductive stage. In contrast, genotype release year was positively associated with maximum net CO2 assimilation rate (Amax), maximum carboxylation rate of Rubisco (Vcmax), maximum electron transport rate (Jmax), seed number, and seed weight, while Amax and Vcmax were positively correlated with seed number and seed weight. These findings indicate that the greater photosynthetic capacity of modern genotypes was not accompanied by faster photoprotective response. Thus, photoprotective regulation has not kept pace with gains in photosynthetic capacity under field conditions. We conclude that rapid NPQ relaxation remains an important target for synchronizing photoprotection with the high photosynthetic capacity of modern soybean lines.

Ecologists remotely track movement steps of animals (e.g., via global positioning systems) and use step selection functions to study the effect of environmental factors upon their movement decisions. Constructing such functions requires pairing each observed step with some number of unobserved but feasible comparison steps. Larger numbers of comparison steps generally yield better estimates but also incur potentially challenging computational demands. Thus, it is important to determine an appropriate number of comparison steps. No established guidance exists for this decision. Here, we use simulated tracks to assess how many comparison steps are needed, fitting each set of steps to a conditional logistic regression model. We monitor errors in estimated effects for several classes of tracks, identifying the number of comparison steps for which mean relative absolute error in estimated effects is consistently low. By this criterion, 32 comparison steps per observed step are needed for our primary class of simulated tracks. Tracks in more homogeneous landscapes, tracks with shorter mean step lengths, or shorter tracks generally require more comparison steps (ranging from 64 to 128 per observed step) to achieve the same level of accuracy. Longer tracks generally require fewer comparison steps (16 per observed step). These results clearly demonstrate that the number of comparison steps influences how well step selection functions estimate covariate effects and provides initial direction in a research area that currently lacks quantitative guidance. Movement ecologists should take care when selecting the number of comparison steps paired with each observed step because those decisions matter.

The southern sea otter (Enhydra lutris nereis) population at San Nicolas Island, California, has been monitored annually since the translocation of 140 sea otters to the island was completed in 1990. Monitoring efforts have varied in frequency and method across years. In 2017, in accordance with the National Defense Authorization Act for Fiscal Year 2016, the U.S. Navy and the U.S. Fish and Wildlife Service formally initiated a sea otter monitoring and research plan to determine the effects of military readiness activities on the growth or decline of the southern sea otter population at San Nicolas Island. The monitoring program, at its basic level, includes quarterly seasonal surveys of population abundance, distribution, and foraging activity. This report presents data from the program with a focus on the recent three years from winter 2023 through winter (February) 2026. From 2023 to 2026, we measured an 8.1-percent per annum decrease in population abundance (95-percent confidence interval =1.1-14.6 percent), with 106 total individuals counted as of February 2026. Historically, sea otter habitat usage at San Nicolas Island was concentrated on the west end of the island. Between 2017 and 2019, we observed increased seasonal usage of the north and south sides of the island, and in 2020-2022, a large (approximately 30-40 individuals) group of sea otters (raft) took up residence off the east end. During 2023-2026 the east end raft disappeared, and sea otters returned to their historical habitat usage patterns at the west end of the island. Foraging data were collected from summer 2023 to winter 2026 on a total of 461 foraging dives in 32 foraging bouts, and the majority of identified prey on successful dives (n=325) were sea urchins (124) followed by snails (48), bivalves (41) and crabs (23). One lobster and one octopus were also identified among the sea otter prey items. We combined these data with data from 2020-2022 to estimate overall energy intake rates that averaged 7.7 kilocalories per minute (95-percent credible interval =6.6-9.1 kilocalories per minute). These results can be useful to the planning of future monitoring and research of sea otters at San Nicolas Island.

Perennial grains represent a promising frontier for sustainable agriculture, but breeding progress is constrained by the accessibility of genotyping and the difficulty of evaluating complex traits expressed for multiple years after establishment across heterogeneous environments. Phenomic selection may help address these challenges by using inexpensive, scalable, high-dimensional phenotypes collected early in development, although the robustness of such predictions across breeding cycles remains uncertain. Here, we compared genomic selection and phenomic selection across two breeding cycles of Thinopyrum intermedium (intermediate wheatgrass; IWG; Kernza(R)), comprising approximately 2,280 individuals from maternal half-sib families evaluated across multiple field sites and years. We constructed relationship matrices from genomic markers and early-life stage phenomic data, including seed and leaf color (HSV), CropReporter multispectral reflectance and indices, and cycle-specific hyperspectral reflectance sensors. Genomic models provided the strongest predictions on average across all field traits in both cycles. Among phenomic predictors, leaf HSV was consistently the most informative, whereas CropReporter and hyperspectral data showed lower and more trait-dependent performance and seed HSV provided little predictive value. Genomic, leaf HSV, and CropReporter models transferred across breeding cycles with little apparent loss of predictive ability relative to within-cycle validation, demonstrating that their predictive signals were not restricted to a single breeding cycle. Early-life stage leaf HSV emerged as a practical, accessible tool for germplasm thinning and early-stage prioritization in perennial breeding programs. Despite limited similarity among relationship matrices, multi-relationship-matrix models rarely improved prediction beyond the stronger constituent single-relationship-matrix model. Together, these results show that early-life stage phenomic data provide reproducible information about agronomic performance expressed years later, but that predictor complexity and data integration do not guarantee improved prediction.

Antibiotic-induced disruption of plant-associated microbiomes has the potential to alter host health beyond the directly exposed generation, yet whether the effects of dysbiosis are transmitted through the seed microbiome remains unknown. Here, we investigated the intergenerational impacts of streptomycin-induced dysbiosis in tomato (Solanum lycopersicum), demonstrated that seed microbiome transfer (SMT) restores progeny microbiome function and disease resistance, and characterized the underlying physiological and genetic mechanisms. Parental streptomycin exposure altered the composition of progeny rhizosphere bacterial communities, reduced expression of defense-associated genes, and increased susceptibility to Xanthomonas perforans. Suppression of immune gene expression was strongly associated with increased disease severity, indicating that parental dysbiosis impaired progeny plants ability to mount effective immune responses. Transfer of the seed microbiome from healthy plant donors partially restored rhizosphere community composition, reduced disease severity and recovered defense gene expression of three genes. Together, our findings demonstrated that antibiotic exposure microbiome disturbance generates intergenerational legacy effects that influence plant immunity and disease susceptibility and seed microbiome transfer can counteract this dysbiosis across generations.

Introduction: Polyendocrine metabolic ovarian syndrome (PMOS), formerly known as polycystic ovary syndrome (PCOS), is a common endocrine disorder affecting women of reproductive age. Besides reproductive and metabolic disturbances, PMOS negatively impacts psychological well-being and quality of life. Despite available treatment options, there remains a need for safe and effective therapies that improve both clinical symptoms and fertility outcomes. Aim: To compare the efficacy of VAMHA and MYRHA tablet combination therapy with standard non-hormonal therapy in restoring regular menstruation. Secondary objectives included assessment of ovulation, menstrual symptoms, polycystic ovarian morphology, hormonal and metabolic parameters, anthropometric measures, and skin manifestations. Study Design: Open-label, randomized, multicentre, prospective comparative clinical study. Methods: Seventy-one women with PMOS were randomized to Group A (n=37) or Group B (n=34). Group A received VAMHA and MYRHA tablets (2 tablets each), while Group B received Metformin 500 mg plus Myoinositol 600 mg (1 tablet), twice daily for 180 days. Data were recorded in Case Report Forms. Statistical Analysis: Continuous variables were summarized using mean and standard deviation, while categorical variables were expressed as frequencies and percentages. Appropriate statistical tests, including Chi-square, were used. A p-value [&le;]0.05 was considered significant. Results: Significantly more participants in Group A achieved regular menstrual cycles than Group B (31 vs. 22; p<0.05). Ovulation occurred in 16 participants in Group A compared with 6 in Group B (p<0.05). Both groups showed significant improvement in menstrual irregularity and related symptoms. Significant reductions in Anti-Mullerian Hormone (AMH), fasting insulin, and body mass index (BMI) were observed in both groups (p<0.05). Resolution of polycystic ovarian morphology occurred in 13 participants (38.23%) in Group A and 10 (33.33%) in Group B. Both treatments were well tolerated with no major safety concerns. Conclusions: VAMHA and MYRHA combination therapy was superior to standard non-hormonal therapy in improving menstrual regularity and ovulation. It also produced favourable metabolic, hormonal, and ultrasonographic outcomes, suggesting its potential as a safe and effective option for comprehensive PMOS management and fertility enhancement.